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1 Introduction 



In this paper we consider the forward and inverse scattering problems on the 
line for operators 

Lfc = + A;V - 9, D = d/dx. (1.1) 

For p = 1, (|1 . 1|) is the Schrodinger operator and k is the wave number, or 
momentum. As is well-known, the scattering data of ( |1 . 1| ) evolves linearly as 
q evolves according to the nonlinear Korteweg-deVries flows. 

For q = and p ^ 1 as a; — ±oo ( |1 . 1| ) is the acoustic scattering problem. 
In this interpretation, k is the frequency. The associated hierarchy of flows, 
of which the first is 

(A = >f-), (1-2) 



2 Vp. 

were introduced by Martin Kruskal WM and attributed to an unpublished 



paper of Harry Dym. 

For g = 1/4 and p — > at infinity, ( |1 . 1| ) was introduced by Camassa 
and Holm [Q, 0, in connection with a nonlinear shallow water model. The 
Camassa-Holm equation itsellQ, in the normalization determined by ( |1 . Ij ) 
with g = 1, is 

1 3 11 

(1 - -D^)ut = -iu^)x - g(MDx - -iuu^x)x, (1-3) 

/ = 2(1 - \d^)u. 

For this equation, considered on the line, it is natural to assume that the 
potential p — > as |x| — >■ oo; the scattering problem on the line is therefore 
singular. The case of periodic potentials has been resolved by Constantin and 
McKean . Camassa and Holm , and Camassa, Holm and Hyman have 
also obtained Blow-up results for certain initial data and have constructed 
two-soliton solutions by direct methods - that is, methods not based on the 
inverse scattering technique. Constantin and Escher have proven global 
existence theorems for the Camassa Holm equation on the real line for a large 
class of initial data. 



^The equation itself had been obtained previously by Fuchsteiner |Q] by the method 
of recursion operators; but that method does not give the isospectral operator. 
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The scattering problem for when p = 1 is by now well studied and 
completely understood; the inverse problem can be solved by the integral 



equation method of Gel'fand and Levitan |]T0|. In the acoustic problem, 
p — ^ 1 as |a;| — >• cxD and g = 0, it is well-known that the classical Liou- 
ville transformation takes Lk to the Schrodinger operator Q. The inverse 
transformation from the Schrodinger to the acoustic problem requires the in- 
version of a differential equation. We show in this paper that this step can be 
reduced to a simple linear quadrature. This greatly simplifies the inversion 
problem. 

It is the acoustic problem, rather than the Schrodinger equation, that is 
relevant to physical applications in elasticity and reflection seismology 



Although the scattering problem for ( |1 . 1|) when p — >■ at infinity and 
g = 1 is a priori singular, a particular change of variables reduces it to the 
Dirichlet problem 

^ + k^pliO] A;) = 0, -1 < e < 1; ^(±1, k) = 0. 



This is a density problem on a finite interval |T2[, |T3|. A second Liouville 
transformation converts this to a Schrodinger problem on a finite interval. 
Thus the inversion problem reduces to solving an inverse Sturm-Liouville 
problem on a finite interval [|ll|, |jl5|. 



2 Acoustic Scattering 

The equation 

(D^ + A:V)^/' = 0, -oo<a;<oo, (2.1) 
arises in scattering problems for the wave equation 

utt - (?{x)u^x = 0, 

where the sound speed c is 1/p. Therefore it is natural to assume p to be 
bounded away from and to have a finite limit as |a;| — > oo. We call ( |2.1| ) 

^We thank Fadil Santosa, of the University of Minnesota, for his helpful information 
regarding the applied literature in acoustic scattering problems. 
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the acoustic scattering problem. We assume throughout that p is real and 
positive, while p — 1 belongs to iS, the Schwartz class of rapidly decaying 
functions. 

We begin by constructing the wave functions of normalized at ±00, 
by the WKB method. We write the wave functions in the form 

ip+{x, k) = i+{x, k)e-'^^^''\ k) = m+{x, k)e'^^^^''\ 

where and m+ are normalized to be 1 at —00 and +00 respectively. Sub- 
stituting this form of the wave functions into equation we find that m+ 
satisfies the differential equation 

m" + 2ikS'm' + [ikS" + k'^{p^ - S'^)]m = 0. 

The term involving k'^ is eliminated by requiring S{x) to satisfy the eiconal 
equation 



We take 



Then 



where 



S'{x)=p{x), S{x) = x+ [p{y)-l]dy. 



S{x) 



X + 0(1) X —00 
a; + 7 + 0(1) X ^ 00 



1=1 [piy) - 1] dy. 

J —00 

We now have 

m" + 2ikpm' + ikp'm = 0, m — > 1 as x — > 00. (2-2) 

The solution to this differential equation is constructed by converting it to 
the Volterra integral equation 

mfx, fc) =1 + / ; (2ik(l — p)m' — ikp'm) dy 

L 2ik 



G{x,y,k)m{y,k) dy, (2.3) 
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where 

G{x,y,k) = 2ike^^^^y-^\p{y) - 1) + ]^{e^^^^y-) - l)p\y). 

This integral equation can be solved for any k in the upper half plane '^sk > 
by the method of successive approximations, since p — 1 and p' belong to 
£^(]R). The solution is analytic with respect to k; we denote it by m^{x, k). 
Wc denote the Schwarz reflection of m+ to the lower half plane by m_(a;, k) = 
m+(a;, k), and extend ip+ to the lower half plane accordingly: 

iP4x, k) = ij+{x, k) = m_(a;, k)e-'^^^''\ 

Similarly the function ^_ defined by 

A;) = 1 - r G{x, y, k)i4y, k) dy. (2.4) 

J —oo 

is analytic with respect to k in the lower half plane and has a Schwarz reflec- 
tion £+{x, k) analytic in the upper half plane. We denote the corresponding 
wave functions by by 

^p±{x,k) ^ £±{x,k)e^'^^^''\ 
The boundary values on the real line satisfy the asymptotic conditions 
lim e^^'^^^^Vi = lim e^^'^^^^Vi = 1- (2-5) 

x-*—oo x— »oo 

Note that 

£+{x,k) =i_{x,-k), m_{x,k) =m+{x,-k), keR. (2.6) 

Moreover G{x, y,0) — so 

4(a;,0) = m±(x,0) = 1. (2.7) 

The scattering data for the acoustic equation is deflned just as it is for 
the Schrodinger equation. For real non-zero k each pair of wave functions 
ip±, Lp± is (generically) independent and thus constitutes a fundamental set 
of solutions. Therefore 

(p+{x, k) = a{k)ilj_{x, k) + b{k)ilj+{x, k) (2.8) 
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for some functions a and b. For real k we use ( [2 .61 ), the definitions of the 
wave functions, and the hmits, to see that (12. 81) becomes 



i+{x, k) ~ a{k) + 6(A;)e2^'=^(^) ~ a{k) + h{k)e^'^^''+^^ as x +oo. 
In view of ( |2.4| ), ( ^.61 ) and the form of G, we obtain for real k that 

1 

a{k) = l + - p'iy)£4y,k)dy; (2.9) 

/OO 
g-W,)|2,fc(p(2/) _ 1) _ ip'(j/)}4(2/, A;)ciy. 
'OO 

It follows from (U) and that 

a(0) = 1, 6(0) = 0. (2.10) 



The formula for a{k) in ( p. 9]) extends to the upper half plane. We take 
the Wronskian of both sides of (|2.8|) with iIj± and obtain 



In fact, taking asymptotics as a; ^ cxd, 

W{iIj_,-^+) = W{e-'''^,e'''^) = 2ikS' ~ 2tk = 2ikp ~ 2ik. 
The expression ( p.ll| ) for a{k) also extends to the upper half plane. 

Theorem 2.1 The reduced wave functions m^^ are analytic in the upper 
half plane, andm_, i_ are analytic in the lower half plane. 

The function a{k) is analytic in the upper half plane. Moreover a(0) = 1 
and a has no zeros. 

Proof: We have proved everything but the last statement. As in the case 
of the Schrodinger equation, the acoustic scattering data for real k satisfy 



a{k) = a{-k), h{k) = h{-k)- |a(A;)|'- |6(A;)|' = 1. (2.12) 

Therefore a has no real zeros. According to (|2.11| ) a zero at k — xuj in 
the upper half plane would correspond precisely to a bound state: an LP' 
wave function with negative energy. This would be an eigenfunction for the 
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operator L = p ^D^ with eigenvalue —k"^ = u"^ > 0. However L is selfadjoint 
and negative in L^(]R, (ix), so it cannot have such an eigenvalue. q.e.d. 

We define the refiection coefficient r by 

rik) = ^ ^k = 0. 

In the absence of bound states, r constitutes the complete scattering data 
for the problem. In fact from ( ^^.121 ) we have 



1 — |rP 

Since a is analytic in Q k 0, tends to 1 as /c — >■ oo, and has no zeros, arg a 
can be recovered from log |a| on the real axis by the Hilbert transform: 

argaW = ip.V. f 

TT J_^ k - t 

Then log a is obtained for ^ k ^ hj the Cauchy integral representation. 



3 The Liouville Transformation 

The acoustic equation may be transformed to the Schrodinger equation by 
the well-known Liouville transformation. By our assumptions, S{x) is a 
monotone increasing function on the line, hence we may define a change of 
variables by ^ = S{x). The variable C, corresponds physically to the time of 
travel. By the chain rule, 

d d d d^ d d 

dx dx d^ ^ d^^ dx'^ ^ d^^ d^ 

To keep track of the relevant variables we define ps{C) by 

Ps(0=P(a;), at^ = 5'(x). 

The mapping / f oS'^ is a unitary map from L^(M, dx) to L^(M, ps d^). 
It carries the negative selfadjoint operator p^'^D'^ to the operator 
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To complete the Liouville transformation we use the unitary map / f yfpl 
from L^(M, psci^) to L^(R, ci^). The corresponding gauge transformation 
(conjugation by p^"^/^) takes the preceding operator to the Schrodinger op- 
erator 

^|-g(0, (3.1) 

where 

Since = p~^Dx = p~^D, we can also express the potential as 

'(^(^))4f-!(^)^^M-^-'^(7p))- 

We remark that 

q{Six)) = ^{^,x} 
where x} denotes the Schwarzian derivative 

D-'^ 3 fD'^ 



2 

We have seen that the transformation from the Schrodinger potential q 
to the associated acoustic potential p is given by the equation (^]2|), together 
with a change of variables. We show below that this transformation may be 
computed by a linear operation on scattering data and a simple quadrature. 



Theorem 3.1 The normalized wave functions (ps,±, 4's,± for the operator 
(|3.1| ), (|3.2| ) are related to the wave functions for the acoustic operator by 

^sAO = ^/p¥)MS-'0; i^s,± = v^MS-'O- (3.4) 

The scattering data (a, b) of the Schrodinger problem is precisely that obtained 
for the acoustic problem ( |2.1| ). 
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Proof: The relation in ( ^.4| ) simply implements the two unitary transforma- 
tions, change of variables and change of gauge. Therefore functions (ps,±, ips,± 
are wave functions for the operator ( |3.1| ). Moreover they have the correct 
asymptotics and analyticity properties. q.e.d. 



Theorem 3.2 The image of the acoustic problem under the Liouville trans- 
formation consists of all Schrddinger operators with Schwartz class potentials 
and no bound states, such that a(0) = 1. 

The inversion of the Liouville transformation is given by 

p,(o = n^, o)^ ^ = ^ + / {jj^ - (3.5) 

Proof: Theorems |2.1| and lOI imply that operators in the range of the Liouville 



transformation have no bound states and have a(0) = 1. The transformation 



equation (|3.2|) shows that the potential is of Schwartz class. The acoustic 
potential p may be recovered from the Schrodinger potential q as follows. By 
(|3.4| ) the wave functions at = are related by 

On the other hand we have observed that the normalized acoustic wave func- 
tions at k = are identically 1. It follows that 



Therefore to reconstruct the function x = x(^) from the Schrodinger poten- 
tial, we can compute ipsi-jO) and use 

dx 1 1 



Suppose, conversely, that g is a Schwartz class Schrodinger potential with 
no bound states, such that a(0) = 1. Then ■?/'s(-,0) = ips{-,0) is real and 
asymptotically 1 in each direction. We prove that iIjs{C,,0) = m'*(^,0) has no 
zeros when there are no bound states. First, note that m^_^_{-,iuj) is real and 
converges uniformly to m^(-,0) as u decreases to 0. If the latter function 
had any zeros they would necessarily be simple, since is a solution of a 
second order differential equation, and therefore mi(-,2u;) would have zeros 
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for small lo > 0. This function converges uniformly to 1 as e +00 and the 
zeros remain simple, so there would be a value u > for which has limit 
as a; — *■ 00. But the corresponding ip^ would be a bound state. Therefore 
'?/'s(-,0) is positive. Because of this we can use ( |3.5| ) to construct a change 



of variables and a potential p. The assumption that q belongs to S implies 
that ips{') — 1 belongs to S. Therefore p — 1 belongs to S. Clearly — g is 
the Liouville transform of the acoustic operator associated to p. q.e.d. 

We have shown that the density p in the time of travel coordinate ^ is 
obtained immediately from the Schrodinger wave function at k = 0. The 
latter can be obtained directly from the Gel'fand Levitan kernel K{x,y), 
since the wave functions are given by 

POO 



4 The Harry Dym Flows 

The Harry Dym flows are related to the operator L = p~'^D'^ exactly as the 
KdV flows are to D| — q. They are derived from commutator conditions 
[L, A] = where A is the Liouville transformation of an operator that deter- 
mines one of the KdV flows. The computation of A is rather complicated, 
and seems to shed no light on the issue. It is therefore more convenient to 
work with Lk = p^{L + k"^) instead. Then the commutator condition must 
be modified to 



[Lk,A] = BLk (4.1) 
where Lk, A, B act on functions ip = ipix, t, k), and Lk and A are given by: 

j=0 

The coefficients fj are taken to be polynomials in k'^. The significance of 
(|4.1| ) is that [L^, A] vanishes on the wave functions, i.e. solutions of L^ip = 0. 

•^By the same argument one can prove that the number of zeros of m*(^, 0) is precisely 
equal to the number of bound states. If there are n bound states, then by oscillation 
theory m'*(^, iujn) has n zeros. The number of zeros of to* is constant as uj decreases from 
ojn to zero; so to*(^,0) also has n zeros. 
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The conditions relating the coefficients of L^, and A may be obtained from 
cross-differentiation of the pair of equations Lkip = 0, Aip = 0. 

Any operator product of the form CL may be added to A without affecting 
( |4.1| ). Therefore even powers of D may be ehminated in favor of even powers 
of k, and it is enough to take n = 1: A = —dt + fD + g. Then 

[L, A] =2f^D^ + + 2g^)D + [k^p^ - ef{p\] + g.. 

+ 2go)D + [k\p\ - ef{p'), + g,, - 2k^p'f,] + 2f,D L. 
The coefficient of D must vanish, so we take 

g = -^Jx- 



Then the conditions for ([4.1|) become 



Setting 



k\p')t = \h.. + k%p')xf + 2p'f,] 

= ^ fxxx + '2p{pf)x- 



i=i 



substituting this expession into the previous identity, and comparing coeffi- 
cients of powers of k, we find that F„ is a constant multiple of 1/p, Fq = 0, 
and the remaining coefficients can be determined from the recursion relation 

4pD(pF,_i) = -D^F,. 

The flow of p is 



For n = 1, we take Fi = 1/p and obtain ( p. .21) . 

Theorem 4.1 Under the Harry Dym flow, the scattering data a{k,t) and 
b{k, t) evolve according to 

d = 0, b = 2ik%. (4.2) 
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Proof: The commutator relation ( |4.1| ) implies that the kernel of L is invariant 
under A; that is, 

Lip{x,t,k)=0 =^ LAip{x,t, k) = 0. 
The wave function 

ip+{x, t, k) = i+{x, t, fc)e-^'=^(^) ~ e-*^^ X -oo. 
On the other hand, for the Dym equation itself, 

A = -dt + —D --D\-\ dt + k'^D, X ±00. 



P 2 Vp. 
Therefore, 

A^+{x, t, k) ~ {-dt + k'^D)e-'^'' = -iPe'^'^^ x -00. 

Since the wave functions are uniquely determined by their asymptotic be- 
havior as a; —i> —00, we conclude that 

A(p+ = —ik^Lf^. 

Similarly, we find that 

Aip± = ±ik^ip±, A(p_ = ik^ip_. 

Therefore, applying A to ( |2.^ ), we obtain 

Aip^ = — ik^Lp^ = —ik^{aip_ + lnp+) 
= — aip- — hip+ + aAip- + bAxfjj^ 
= — dip- — bip+ — ik^atp- + ik^b'ip+; 



and equations ( ^4.2|) follow immediately from the independence of ip+ and ip-. 
q.e.d. 

5 The Camassa-Holm spectral problem 

We consider the operators 

L = D^ + ep^-l, (5.1) 
12 



where D = d/dx, introduced by Camassa and Holm We are interested 
in a singular case of (|5.1|): p is positive but p ^ rapidly at infinity. As a 
consequence, the natural normalization for the wave functions is independent 
oft. 

lim e"^v9(x) = 1 = lim e'^tjj{x). (5.2) 
The Camassa-Holm equation (|1.3|) implies that the evolution of is given 



by 

(p^)t = u(p^)^. + 2M^p2. 
Under our assumption that p is strictly positive, we have the equivalent form 

pt = (up),, (4 - D^)u = 2p\ (5.3) 

Let us assume that 

f oo 

e^\y\ p{yf dy <oo. (5.4) 

) 

Then, from the second equation in ( |5.3| ), 

1 

u{x) = - e-^\--y\p{yfdy. 

^ J -oo 

It follows that 

/oo 
e^^yp{yfdy. (5.5) 
-CXD 

Therefore the evolution ( ^.3|) is consistent with the assumption (5^), and 
with the stronger assumption p(x) = 0(e~^'^') as |x| oo. 

For reasons that will become clear, a natural class of functions p to con- 
sider are positive C°° functions that have asymptotic expansions 

oo 

p(x) ~ a^e"^'^'^' as X ^ ±oo, a]" = aj,*" > 0. (5.6) 

It follows that u has similar expansions. The leading terms in these 
expansions are given by ( |5.5| ) and the remaining terms can be computed 
from ( ^.6|) and the second equation in (|5.3|) . These expansions are consistent 
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with the evolution ( |0| ) and the time dependence of the coefficients can be 
computed from these equations. 

We observe next that for a particular choice of potential p the Liouville 
transformation trivializes the operator p^'^{D^ ~ !)• By ( p.3| ), the latter 
transforms to a Schrodinger operator with potential 



q{S{x)) 



ID^p _ 3/Dp\ '- 



+ 1 



Therefore if we choose 



1 



Pq{x) 



cosh X 



sech^x 



(5.7) 



the potential go vanishes. Note that Pq{x) has an expansion ( ^.6|) . Since 
sech^x = D tanh x, we can take as transformed variable C, = So{x) = tanh x. 
Then the range of the transformation is the interval — 1 < < 1.^ 

The positive smooth potentials p that have asymptotic expansions 
are precisely the potentials that can be written in the form 



p(x) = (yf(tanh x) po{x) = 5f(tanh x)sech^x, 
geC^i[-l,l]), g>0, g{-l) = gil). 



(5.8) 

(5.9) 



The Liouville transformation generated by po{x) brings Po^(Z)^ + k'^p'^ — 1) 
to the form 



(5.10) 



This is the density problem for a string on the finite interval |jT2[, 
The corresponding wave functions are related by 



'4's{C,k) = \/pJX)'ipa{x,k) = ^/p{x)^aix,k) = sech x^a{x , k) . (5.11) 



^It was already remarked by Camassa and Holm ||^ that the spectrum of the operator 
(5.1) is discrete when p —^ sufficiently rapidly as x ±oo; and they analyzed the 
special case when Pq — ylsech^(a;). Note, however, that it is the potential p§ — sech''(x) 



that trivializes under the Liouville transformation. 
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Hence the corresponding eigenf unctions vanish as ( — > ±1, and the appro- 
priate boundary conditions for ( |5.10|) are Dirichlet boundary conditions. 

The Liouville transformation from with d^/dx = p and the Liouville 
transformation from ( p.lOl ) with d^/d( = g{() arrive at the same point: a 
Schrodinger operator on a finite intervaL The length of the interval is 2M 
where 

M = lj^ p{x) dx = ^ j'^ g{0 dC (5.12) 

and we normalize the interval to be (— M, M). The expression for the poten- 
tial g(^) in terms of g is given by ( p.3| ): 

<iim) = i^,{ix}. (5.13) 



6 Camassa-Holm scattering data 

In this section we describe the scattering data for ( p.l| ) and find its evolution 
under the Camassa-Holm flows. 

We have seen that the Liouville transformation is a unitary equivalence 
between the negative operator p~'^{D^ — 1) and the Schrodinger problem on 
a finite interval. Under this transformation the normalized wave functions 
( |5.2| ) are multiplied by p, which vanishes at oo. On the finite interval, the 
asymptotic conditions become Dirichlet conditions: vanishing at an endpoint. 
Thus it is natural to take as scattering data for ( [5.1| ) the Dirichlet spectrum 
and the associated coupling coefficients. Note that the wave functions for 
( |5.1CI| ) at = have the form a( + 6; the corresponding wave functions 
for the Schrodinger operator do not satisfy Dirichlet conditions. Thus the 
Dirichlet eigenvalues are strictly negative. 

At the eigenvalues —k"^ the wave functions (p and ip are linearly dependent: 
there is a coupling coefficient such that 

Thus we take as scattering data for p the countable set: 

{kn,Cn}. (6.1) 
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The Camassa-Holm equation is obtained from the commutator relation- 
ship (|4.1|), with L given by ( ^.1|) and 

A=-— + aD--a^, a = u{x,t) + —. 

Note that the interval (— M, M) remains constant, in view of (|5.12| ) and ( ^.3|) . 

To determine the evolution of the coupling coefficients, we proceed as in 
the Harry Dym flow, with suitable modifications. The compatibility relation 
(|4.1| ) implies that the kernel of L is invariant under A; hence Aip{x, t, k) is a 
linear combination of y^{x,t, k) and ip{x,t, k). The exact linear combination 
is determined by evaluating the asymptotic behavior of Aip as x — > ±00. 
Under our assumptions, u and tend to zero exponentially as x tends to 
—00. Again, the leading asymptotics of ip are independent of t, so we find 
that the leading asymptotics of Aip as x — — 00 are 

1 Dc^ 6^ 

A(p = {—dt + aD cLx)(p 



hence 



Similarly, 



At a bound state, 



2 ^'"^ F f.2 



Aif{x,t,k) = ^!f{x,t,k). 



Aip{x,t,k) = -^ij{x,t,k). 



ip{x, t, kn) = Cn{t)^{x, t, kn), 



hence 



i'n = fn = Aifn = A{Cn1pn) 



-Cnlpn + CnAlpn) = { -C„ - ) Vn- 



Since ipn is not identically zero, we have 

^"(^) = -IX' ^n{t) = e-2*/^"c„(0). (6.2) 

We next turn to the characterization of the image of p^^{D^ — 1) under 
the Liouville transformation to a Schrodinger operator. 
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Theorem 6.1 The necessary and sufficient conditions that the Schrddinger 
operator D| — q on the interval (— M, M) be in the range of the Liouville 
transformation associated to ( p.l| ), where p = g{tanh x)secf?x and g is a 
positive smooth function on [—1, 1] with (?(— 1) = giX), are that q he smooth 
on [— M, M] and that the Dirichlet spectrum of — q be strictly negative. 

Proof: It is clear from the formula that describes q in terms of go that q 
must be smooth, and we have noted above that the Dirichlet spectrum must 
be strictly negative. 

Conversely, suppose that the Schrodinger potential q is smooth and the 
Dirichlet spectrum is strictly negative. By standard oscillation theory the 
non-zero solutions to 

iDl-q)ipo = 0, M-M)=0, (D|-g)^o = 0, MM) = 

do not change sign on the interval. Therefore the (unique) solution to the 
Dirichlet problem 

{Dj - q)ip = 0, -^{-M) = iP{M) = 1 

is a linear combination of ipQ and ipQ that is strictly positive on the closed in- 
terval. The inverse Liouville transformation must be associated to the change 
of variables dQ/dS, = X/ipl^)"^, where A is determined by the normalization 
condition 
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